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Why a falling drop does not in genera 
behave like a rising bubble 

Manoj Kumar Tri path i, Kirti Chandra Sahu & Rama Govindarajan* 

Department of Chemical Engineering, Indian Institute of Technology Hyderabad, Yeddumailaram 502 205, Andhra Pradesh, India. 

Is a settling drop equivalent to a rising bubble? The answer is known to be in general a no, but we show that 
when the density of the drop is less than 1.2 times that of the surrounding fluid, an equivalent bubble can be 
designed for small inertia and large surface tension. Hadamard's exact solution is shown to be better for this 
than making the Boussinesq approximation. Scaling relationships and numerical simulations show a 
bubble-drop equivalence for moderate inertia and surface tension, so long as the density ratio of the drop to 
its surroundings is close to unity. When this ratio is far from unity, the drop and the bubble are very 
different. We show that this is due to the tendency for vorticity to be concentrated in the lighter fluid, i.e. 
within the bubble but outside the drop. As the Galilei and Bond numbers are increased, a bubble displays 
underdamped shape oscillations, whereas beyond critical values of these numbers, over- damped behavior 
resulting in break-up takes place. The different circulation patterns result in thin and cup-like drops but 
bubbles thick at their base. These shapes are then prone to break-up at the sides and centre, respectively. 

Ask a fluid mechanician the following question: If a bubble of density pi rises in a fluid of density p 2 > would a 
drop of density p 2 rise in fluid of density p x in the same way, if gravity were reversed? The answer is likely 
to be a yes (except from experts on the topic), with a caveat that viscosity must be appropriately reassigned. 
The correct answer is a no 1 , and in this paper we study the physics behind this difference. We show how vorticity 
tends to be concentrated in the lighter fluid, and this changes the shape and the entire behavior. By the term 
bubble we mean a blob of fluid which is lighter than its surroundings, and by drop a fluid blob that is heavier than 
its surroundings. We discuss the limit in which a drop and bubble can behave the same, and present two 
alternatives in this limit to design an equivalent drop for a bubble. 

Due to their importance in industry and in natural phenomena, and due to their appeal to our scientific 
curiosity, bubbles and drops have been studied for a long time. The literature on this subject is therefore huge, with 
several review articles by several authors 2 " 4 and others. Appealing introductions to the complexity associated with 
bubble and drop phenomena can also be found in Refs. 5, 6. There are yet many unsolved problems, which are the 
subject of recent research (see e.g. 7 " 13 ). Bubbles and drops have often been studied separately, e.g. 11,14 " 18 on bubbles 
and 19 " 24 on drops. However, there is also a considerable amount of literature which discusses both together, 
e.g. 25-30 . 

When the bubble or drop is tiny, it merely assumes a spherical shape, attains a terminal velocity, and moves up 
or down. An empirical formula for the terminal velocity of small air bubbles was found by Allen 31 in 1900. Later, 
two independent studies by Hadamard 25 and Rybczynski 32 led to the first solution for the terminal velocity and 
pressure inside and outside of a slowly moving fluid sphere in another fluid of different density and viscosity. The 
spherical vortex solution due to Hill 33 has been a keystone for most of the analytical studies on the subject. Later 
studies 1 showed that at low Reynolds and Weber numbers, drops and bubbles of same size behave practically the 
same way as each other, both displaying an oblate ellipsoidal shape. 

Bigger drops and bubbles are different. A comparison of bubble and drop literature will reveal that in the typical 
scenario, bubbles dimple in the centre 14 ' 34 while drops more often attain a cup -like shape 13 ' 35 . This difference 
means that drops and bubbles which break up would do so differently. The dimples of breaking bubbles run deep 
and pinch off at the centre to create a doughnut shaped bubble, which will then further break-up, while drops will 
more often pinch at their extremities. A general tendency of a drop is to flatten into a thin film which is unlike a 
bubble. There could be several other modes of breakup (see 35 ) like shear and bag breakup for drops falling under 
gravity and catastrophic breakup at high speeds. 

We mention a few significant numerical studies before we move on. Han & Tryggvason 35 studied the secondary 
break-up of drops for low density ratios, whereas Tomar et al. 56 investigated the atomization involving the break-up 
of liquid jet into droplets. We also note that recent developments in the field of lattice Boltzmann simulations 54 ' 55 
allow one to simulate multiphase flows with high density ratio by using a clever doubly-attractive pseudo -potential. 
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Figure 1 | Schematic diagram showing the flow domain as a box 
surrounding a bubble/drop and the axis of symmetry as a dot-dashed line. 

The outer and inner fluid phases are designated as V and Y, respectively. 

We now pose the problem and discuss our approach to its 
solution. 

Problem formulation and approach. The geometry and equations. 
The drop or bubble is initially spherical with radius R and has 
viscosity and density and p z - respectively, which are different 
from the respective quantities \i 0 and p 0 for the surrounding fluid. 
The fluids are assumed to be immiscible, incompressible and 
Newtonian. Figure 1 is a schematic diagram of the initial configura- 
tion, where x and z represent the horizontal (radial) and vertical 
(axial) directions, respectively. The subscripts Y and V stand for 
inner and outer fluid, respectively, a is the surface tension 
coefficient between the two fluids and g the acceleration due to 
gravity. We study an axisymmetric flow, with the velocity and 
pressure fields denoted by u(u, w) and p respectively, wherein u 
and w are the velocity components in the x and z directions, 
respectively. We treat the flow as isothermal. 

Given the physical quantities in the problem, the Buckingham-pi 
theorem requires four independent dimensionless parameters to 
completely describe the system and the ones we choose are the den- 
sity ratio r, the viscosity ratio m, the Galilei number Ga and the Bond 
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Figure 2 | Vertical location of the center of gravity as a function of time 
for a drop (r = 1.214 and m = 76), an equivalent bubble based on 
Boussinesq approximation (r = 0.786 and m = 76) and an equivalent 
bubble based on conditions (17) and (18) (r = 0.85 and m = 0.1). The rest 
of the parameters are Ga = 6 and Bo = 0.0005. 
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Figure 3 | Evolution of (a) bubble shape with time for r = 0.9, mix = 0.56. 
(b) drop shape with time for r = 1.125, mix = 0.56. (c) drop shape with 
time for r = 1.125, with viscosity obtained from equation (17). The 
direction of gravity has been inverted for the drop in order to compare the 
respective shapes with those of the bubble. In all three simulations, Ga = 
50, Bo = 50, and the initial shape was spherical. 

number Bo. These are defined respectively as r = p,-/p 0 , m = \ij\i 0 , Ga 
= p 0 g V2 R 3/2 l p 0 and Bo = p 0 gR 2 /a. The governing equations are thus 
non-dimensionalised using R, \/gR, fi 0 an d p 0 as length, velocity, 
viscosity and density scales, respectively. Note that the Galilei num- 
ber is like a Reynolds number, but defined by the gravity velocity 
scale rather than the actual forward velocity of the bubble. Since the 
latter is unknown a priori, and depends on several physical quant- 
ities, the present velocity scale is a cleaner one to adopt, and we do so 
in our study. For our discussions however it will be useful also to 
define a Reynolds number Re = p o U 0 R/fi o and a Weber number and 
We = p 0 UqR/(t, where U 0 is a typical bubble velocity. Note that each 
of the four dimensionless parameters, when independently varied 
can change the dynamics of the blob. Along with the continuity 
equation: 



V-u = 0 9 



(1) 



the momentum equation can be written with a continuum model for 
surface tension forcing (see 36 ) as 



Du 



Vp 



Vn 



+ 7F: V '(^(Vu + Vu r )) - j+ — <5(x-x S )n, (2) 



pGa 



pBo 



wherein; is the unit vector in the vertically upward direction; 3 is the 
Dirac-delta function, x s a location on the interface and n its unit 
normal and D[ = d/dt + w-V] is the material derivative. On the right 
hand side, we have the pressure, viscous, gravitational and surface 
tension forces per unit volume respectively. Note that p and /n are 
now dimensionless quantities, respectively scaled by the density and 
viscosity of the outer fluid. 

Thus the surrounding or outer fluid obeys 

DU °--Vp 0 +-^V 2 u o -j, (3) 



Dt 



Ga 



whereas the inner fluid, namely the drop or bubble, satisfies 



Dui 

~Dt 



Vpi 



r rGa 



(4) 



with the interfacial conditions at x = x s being continuity in velocity 
and the relevant stress balance. The surface tension term appears in 
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Figure 4 | Evolution of (a) bubble and (b) drop shapes with time, when 
densities of outer and inner fluid are significantly different. As before, for 
the drop (right column), the direction of gravity has been inverted. In both 
simulations Ga = 50 and Bo = 10. The other parameters for the bubble 
system are r = 0.5263 and m = 0.01, while for the drop system r = 10 and 
m = 0.19. Note the shear breakup of the drop at a later time. Shown in 
color is the residual vorticity (Kolar, 2004). 

the normal stress balance, which ispi\ s = p Q \ s + kIBo wherein k is the 
local surface curvature defined as k = V i 

Numerical simulations. In the framework of the diffuse-interface 
method 37 , we express the surface tension term in the momentum 
equation (2) in the form 

pBo' 



(5) 



and solve it along with the Cahn-Hilliard equation 



(6) 



and the continuity equation (1). Here c is the volume fraction of the 
outer fluid, M = c( 1 — c) is the mobility, ^) = e~ l ooW' (c) — e<7(xV 2 c is 



the chemical potential wherein e is a measure of interface thickness, 
1 2 

*F(c) = - c ( 1 — c) is the bulk energy density and a is a constant; Pe 

= RVI{M c (j) c ) represents the Peclet number wherein M c and <j> c are 
fixed characteristic values of mobility and chemical potential, 
respectively (see 37 and references therein for more details). The 
dimensionless density and viscosity of the fluids are taken to be linear 
in c, as 



p = c + r(l — c), fi = c-\-m{l — c) 



(7) 



where the value of c is one and zero for outer and inner fluids, 
respectively. 

Symmetry boundary conditions are applied at x = 0 and Neumann 
boundary conditions are used at the other boundaries. To avoid 
numerical effects from the outer boundaries on the flow a computa- 
tional domain of width 16.R and height 64R is found to be sufficient 
for most of the cases, consistent with the study of Uno & Kintner 38 . 
Where needed, larger domains have been used, and it has been 
checked that the results are independent of the domain size. In addi- 
tion, the open source code Gerris, created by Stephaen Popinet 39 has 
been used for simulating drops which become very thin as compared 
to their initial radius. The adaptive grid refinement feature in Gerris 
was exploited for our benefit which saved a lot of time in our numer- 
ical study. The numerical method used and its validation are discussed 
in detail in the supplementary material. This code has previously been 
used extensively to study various miscible and immiscible flows (see 
e.g. Mishra et al. 40 ). The supplementary material contains more 
information about the numerical methods and comparisons with 
analytical results in the creeping and potential flow limits. 

Drops versus bubbles: when must we discriminate? Let us first con- 
sider a solid spherical object falling through the fluid. Its equation of 
motion in Stokes flow (small Ga), rendered dimensionless with the 
same scales as before, is 



GaH = 



Ga - —H(\ + 0(Ga)), 



(8) 



where H is the vertical location of the centre of the sphere, and a dot 
represents a time derivative. Maintaining Ga the same, it is easy to see 
that with gravity inverted and setting 



r b = 2-r d , 



(9) 



a light solid bubble would attain the same steady state velocity. The 
subscripts d and b have been added to distinguish drop from bubble. 
Equation (8) makes it obvious however that in the transient phase, 
the accelerations of the bubble and drop must be different, which 
means that at steady state, the vertical locations of the drop and 
bubble will have a constant displacement. Note also that for every 
r b < 1, we can find a realistic r d . However for r d > 2, there is no 
equivalent positive r b . In other words, every solid bubble has an 
equivalent solid drop, but drops heavier than r d = 2 do not have 
an equivalent bubble. 

The next question is, does the same relationship carry across to 
fluid bubbles and drops. The fact that the gradients in the viscous 







Figure 5 | Streamlines in the vicinity of a bubble for t = 1, 2, 3 and 4 for parameter values Ga = 50, Bo = 10, r = 0.5263 and m = 0.01. The shape of the 
bubble is plotted in red. 
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Figure 6 | Streamlines in the vicinity of a drop for t = 1, 2, 3 and 4 for parameter values Ga = 50, £o= 10,r= 10 and m = 0.19. The shape of the bubble 
is plotted in red. 



stresses and pressure inside the drop do not disappear any more 
makes the problem more complicated but, guided by Taylor and 
Acrivos 1 we may look for an equivalence in the limit Re ^ 1 and 
We <C 1, i.e., when we will have small inertia and a near- spherical 
drop. Hadamard's formula 25 for drag will apply in this case, and the 
equation of motion now is 

q(rn 



GaH = 



Ga- 



-H{\ + 0(Ga)). 



10) 



For a small or slowly moving drop, the 0(Ga) terms can be neglected. 
For a Hadamard flow 



q(m) = 



3(2 + 3m) 
2(l + m) : 



(11) 



and it is seen that this expression reduces to the factor 9/2 in equation 
(8) in the limit of infinite inner-fluid viscosity. Equation (10) is a 
first-order ordinary differential equation in H, which can be inte- 
grated analytically. Given that the drop starts from rest, its velocity is 



H d = 



r d Ga d 



-idW t 

l—e r d Ga d 



and equivalently for a bubble with gravity reversed; 
r b Ga b 



q b (m) 





r -%H t ] 


(H 









(12) 



(13) 



For the terminal velocities of the two to be the same, we must have 
r d Ga d ( ^ _ l\ r b Ga b ( 1 
r dj 



q d 



qb 



n 



(14) 



For the outer fluid to have the same dynamics in drop and bubble 
cases, we must have, from equation (3), 



Ga d = Ga b , 



(15) 



while for the pressure jump across the interface to be the same in the 
two cases we need 



Bo d = Bo b . 



(16) 



In addition, we would like the two to have the same history, i.e., to 
have been at the same position at the same time. This requirement is 
met for a Hadamard flow by setting 

r d (l + m d ) r b (l + m b ) 



2 + 3m d 



2 + 3m b 



(17) 



where we have made use of the expression (11). Using equation (17), 
equation (14) yields 

n 



r d -- 



2r b -l 



(18) 



Keeping the Galilei numbers the same, condition (17) will yield 
positive values of m d for a given m b only when r d is less than a certain 
quantity. For the range 0.01 <m b < 100, we find r d must be less than 
— 1.2, and therefore r b >~ 0.86. Thus the equivalence applies only in 



a small range of density ratios. In other words, at low Ga and Bo, a 
drop falling through a surrounding fluid of slightly lower density will 
display dynamics closely matching that of a corresponding bubble 
rising through slightly denser fluid. The dimensionless parameters of 
the drop and bubble are related by equations (15) to (18). On the 
other hand, as seen from equation (18), when r is significantly dif- 
ferent from 1 the dynamics of a drop should be qualitatively different 
from any bubble. 

We now examine whether any equivalence is possible at higher Ga 
and Bo between a drop and bubble. The force balance, written in 
order of magnitude terms (in dimensional form) reads as 



R 



■0 



' R 



+ 0 



i? 2 



+ 0(gAp) + 0 



(19) 



where tildes in the superscript represent dimensional quantities. It is 
clear that for the pressure field in the surrounding fluid in the bubble 
and drop case to be the same, the nonlinear term in the velocity must 
be small, i.e., the flow must be gravity and not inertia dominated, or 
the Froude number, ^/(gR) <C 1. In this case, the drag scales as the 
forward velocity of the drop, so we have V ~ ApgR 2 /in, and we have 
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Figure 7 | Evolution of (a) bubble and (b) drop (gravity reversed) shapes 
with time. Parameters for both bubble and drop systems are: Ga = 100, 
Bo = 50 and m = 10. The density ratio for the bubble and drop are r = 0.52 
and r = 13 respectively, based on equation (18). 
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Figure 8 | Variation of dimple distance versus time for Ga = 50, r = 
7.4734 X 1(T 4 , m = 8.5136 X 10" 6 . 
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-l) 2 Ga 2 <Cl. 



(20) 



Similarly, for surface tension to overwhelm inertial effects as far as 
the dynamics is concerned, and thus not allow for significant shape 
change of the drop, we must have p V 2 /^ <C <j/R 2 , yielding 

1 



(r-l) z Ga 2 < 



rBo 



(21) 



Thus an equivalence between a drop and bubble must be possible 
when r — > 1, even if Ga and Bo are not small, so long as inequalities 
(20) and (21) are satisfied. In this limit, it is usual to apply the 
Boussinesq approximation, by which the fact that the densities of 
the drop and the surrounding fluid are different will only affect the 
gravitational force term in the Navier-Stokes equation. In other 
words, r is set to 1 everywhere except in the gravity term. For drops 
it was noted by Han & Tryggvason 35 that the Boussinesq approxi- 
mation gives reasonable results for 1 < r < 1.6. Interestingly, if we 
had begun by making the Boussinesq approximation, the Navier- 
Stokes equation would yield an equivalent bubble for a given drop 
with r b = 2 — Yd and = m d . We examine in the following section 
which of these choices will work better. 

Results 

Drop and bubble similarities and differences. In all cases where we 
compare drops and bubbles, the Galilei number and Bond number 
are kept the same for the bubble and the drop. In the case where 
inertia is small and surface tension is large, we saw that an equivalent 
bubble for every drop may be found when r is close to 1. We show in 
figure 2 the vertical location of the centre of gravity, 
nondimensionalised with the initial radius of the drop, of small 
inertia and high surface tension, with a density ratio r = 1.214, 
and viscosity ratio m = 76.0 as compared to its surroundings, 
obtained from numerical simulations. Also shown are the vertical 



locations of an equivalent bubble obtained from equations (17) and 
(18), and also one with = 2 — r d , and = m d . For ease of 
comparison, the direction of gravity is reversed for the drop in all 
the results presented. It is seen that the bubble obtained by equating 
the Hadamard drag produces a better equivalent bubble than one 
obtained by a priori making the Boussinesq approximation. All the 
blobs remain practically spherical till the end of the simulation 
shown. 

At small inertia and moderate surface tension (not shown), bub- 
bles and drops do not remain spherical, but with the density ratios 
close to 1, and density and viscosity ratios related by equations (18) 
and (17), a bubble and drop display the same shape and dynamics as 
a function of time. We now consider bubbles and drops of higher 
inertia, where the inequalities (20) and (21) are not obeyed. Drop 
shapes obtained numerically for a density ratio close to unity are 
shown for this case in figure 3(b). The Galilei and Bond numbers 
are the same for the bubble and the drop, and the densities are related 
by equation (18). The Reynolds number based on the terminal velo- 
city of the bubble and the drop is about 16. The kinematic viscosity 
ratio mlr for the bubble is kept the same as the drop in figure 3(a) and 
related by equation (17) in figure 3(c). It is apparent that a drop and 
its equivalent bubble behave qualitatively the same. The velocities of 
the drop and bubble are closer to each other when the viscosity 
relation (17) is used whereas the shapes are closer together when 
they have the same kinematic viscosity ratio. 

When the density ratio is far from unity, no equivalence is pos- 
sible. Equation (17) is no longer valid, nor possible to satisfy. We 
therefore compare drops and bubbles of the same mlr. Figure 4 
makes it evident that neither the shape nor the velocity of the drop 
and bubble are similar to each other. 

Shown in color in this figure is the residual vorticity Q 41 , which is a 
good measure of the rotation in a flow. It is just the vorticity from 
which the shear part has been subtracted in a manner slightly differ- 
ent from that in the Okubo -Weiss parameter, and is discussed in 
greater detail in the supplementary material. Especially at later times 
in this figure, it is evident that residual vorticity is concentrated 
within the bubble but outside the drop. This is the primary difference 
between a bubble and a drop. The region of low pressure and high 
vorticity tends to lie in the lighter fluid. In the case of the bubble, this 
causes an azimuthally oriented circulation in the lower reaches, 
which then leads to a fatter base and aids in a pinch-off at the top 
of the bubble. In the case of a drop, the vorticity being outside means 
that the lower portion of the drop is stretched into a thin cylindrical 
sheet, and an overall bag-like structure is more likely. Also a pinch off 
in this sheet region is indicated rather than a central pinch-off. We 
present in figures 5 and 6 streamlines at various stages of evolution in 
this simulation. Closed streamlines are visible in the region of lower 
density, indicative of regions of maximum vorticity being located in 
the lighter fluid. 

The fact that regions of low pressure and high vorticity would 
concentrate in the less dense fluid follows directly from stability 
arguments. A region of vorticity involves a centrifugal force oriented 
radially outwards, i.e., pressure increases as one moves radially out- 
wards from a vortex. There is a direct analogy (see e.g. 42 ) between 







Figure 9 | Streamlines in and around the bubble at time, t = 1, 1.5, 2.0 and 2.5 respectively, for Ga = 50, Bo = 29, r = 7.4734 X 10 
10~ 6 . The shape of the bubble is plotted in red. 
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Figure 10 | Streamlines in and around the bubble at time, t = 2.5, 5, 7, 9 and 1 1 respectively, for Ga = 50, Bo = 15, r = 7.4734 X 10~ 4 and m = 8.5136 X 
10~ 6 . 



density stratification in the vicinity of a vortex and in a standard 
Rayleigh Benard flow. In the latter, we have a stable stratification 
when density increases downwards. In the former, we have a stable 
stratification when density increases radially outwards, i.e., when the 
vortex is located in the less dense region. 

Figure 7 is a demonstration that for the same outer fluid even if the 
viscosity of the bubble and the drop were kept the same, and only the 
densities of the two were different, the behavior discussed above is 
still displayed. 

Our results thus indicate that density is the dominant factor rather 
than viscosity in determining the shapes of inertial drops and bub- 
bles. In particular, the vorticity maximum tends to migrate to the 
region of lower density, and this has a determining role in the shape 
of the structure. Since large density differences bring about this dif- 
ference, these are effectively non-Boussinesq effects. 

We note that given the large number of parameters in the problem, 
including initial conditions, which we have kept fixed, the location of 
maximum vorticity in the less dense region may not be universally 
observed in all bubble and drop dynamics. For example, the Widnall 
instability 43 in drops resembles the central break-up we have dis- 
cussed. In the usual set-up of the Widnall instability, the densities 
of the inner and outer fluid are close to each other, so we may expect 
the drop to behave similar to a bubble, and the initial conditions are 
not stationary. 

Our arguments above on vorticity migrating to the lighter fluid do 
not depend on gravity being present. We also confirm this in simula- 
tions which obtain the motion of a bubble and a drop started with a 
particular initial velocity in a zero-gravity environment. These are 
shown in the supplementary material. 

Before breakup. Break-up of drops and bubbles is typically a three- 
dimensional phenomenon on which much has been said, see e.g. the 
experimental studies of 44 " 46 theoretical work of 47 ' 48 and numerical 
study of 3 . The transient behavior of liquid drops has been 
discussed extensively 49 " 52 , especially in the context of internal com- 
bustion engines, emulsification, froth-formation and rain drops. 
However, the transient behavior of bubbles has not been 
commented upon as much, and we make a few observations, 
regarding large-scale oscillations, that are not available in the 
literature to our knowledge. We note that since our simulations are 
restricted to axisymmetric break-up they may not always capture the 
correct break-up location or shape. 

Various parameter ranges are covered in numerous papers in the 
past 100 years, and it is known that drops and bubbles break up at 
higher Bond numbers. The Bond number below which a bubble of 
very low density and viscosity ratio does not break, but forms a stable 
spherical cap, is about 8 (see e.g. 53 ) which is similar to that found in 
our simulations. We begin by associating a time scale ratio with the 
Bond number. It may be said that surface tension would act to keep 
the blob together whereas gravity, imparting an inertia to the blob, 
would act to set it asunder. The respective time-scales over which 
each would act may be written as T s = yj pR 3 /<7 for surface tension 
and T g = y/R/g for gravity. The ratio 



T 



1 

Bo 



(22) 



is a measure of the relative dominance. At Bo ^> 1, surface tension is 
ineffective in preventing break-up, and we may expect a break-up at a 
time of 0(1), since we use gravitational scales. For Bo ~ 1, it is 
reasonable to imagine a tug-of-war to be played out between surface 
tension and inertia in terms of shape oscillations, with a frequency of 
0(1). Since a bubble usually breaks up at the centre by creating a 
dimple, the vertical distance D d of the top of the dimple from the top 
of the bubble is a useful measure to observe oscillations. Figure 8 
shows the dimple distance as a function of time for various Bond 
numbers, and our expectations are borne out. 

Figures 9 and 10 are typical streamline patterns in the vicinity of 
bubbles in the break-up and recovery cases respectively. Instan- 
taneous streamlines are plotted by taking the velocity of the foremost 
point of the bubble as the reference, but the picture is qualitatively 
unchanged when the velocity of the centre of gravity of the drop is 
chosen instead. Both cases are characterized by a large azimuthal 
vortex developing within the bubble initially. In the break-up case 
(figure 9), this vortex is sufficient to cause the bubble surface to 
rupture and obtain a topological change, from a spherical-like bubble 
into a toroidal one. In all the cases of bubble recovery we have 
simulated, of which a typical one is shown (figure 10), there develop 
at later times several overlaying regions of closed streamlines, which 
act to counter the effect of rupture by the primary vortex, and to bring 
back the drop to a shape that is thicker at the centre. During each 
oscillation, we see the upper and lower vortices form and disappear 
cyclically. 

A typical breaking drop, with its associated streamlines is shown in 
figure 1 1. As discussed above, the breakup is very different from that 
of a bubble, since the primary vortical action is outside, and causes a 
thinly stretched cylindrical, rather than toroidal shape. 

The vortex in the wake of the drop tends to stretch the interface 
(and surface tension is not high enough to resist the stretching) 
which leads to an almost uniformly elongated backward bag. New 
eddies are formed due to flow separation at the edge to the "bag" and 
a toroidal rim is detached after some time from this bag. A drop too 
responds to Bond number, but the response is shown in terms of an 






Figure 1 1 | Streamlines in and around the drop at time, t = 4.5, 6 and 7.5 
respectively, for Ga = 50, Bo = 5, r = 10 and m = 10. 
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Figure 12 | Variation of break-up time with Bond number for Ga = 50, r 
= 10 and m = 10. 

early break-up at high Bond numbers, as seen in figure 12. The shape 
at break-up too evolves with the Bond number, as shown. 

Conclusions 

A bubble and a drop, starting from rest and moving under gravity in a 
surrounding fluid, cannot in general be designed to behave as one 
another's mirror images (one rising where the other falls). We have 
shown that the underlying mechanism which differentiates the 
dynamics is that the vorticity tends to concentrate in the lighter fluid, 
and this affects the entire dynamics, causing in general a thicker 
bubble and a thinner drop. However, if inertia is small, surface ten- 
sion is large, and a drop is only slightly heavier than its surrounding 
fluid, a suitably chosen bubble can display dynamics similar to it. In 
this limit, the Hadamard solution can be exploited to design a bubble 
with its density and viscosity suitably chosen to yield the same accel- 
eration at any time as a given drop. We are left with an interesting 
result: while a solid 'bubble' can never display a flow history which is 
the same as a solid 'drop', a Hadamard bubble can. Also, although 
density differences are small, the Boussinesq approximation cannot 
lead us to the closest bubble for a given drop. 

We find numerically that a similarity in bubble and drop dynamics 
and shape is displayed up to moderate values of surface tension and 
inertia, so long as the density ratio is close to unity. In axisymmetric 
flow, the vorticity concentrates near the base of the bubble, which 
results in a pinch-off at the centre whereas the cup-like shape dis- 
played by a drop, and the subsequent distortions of this shape due to 
the vorticity in the surrounding fluid, encourage a pinch-off at the 
sides. Bubbles of Ga higher than a critical value for a given Bo will 
break up at an inertial time between 2 and 3. For Ga or Bo just below 
the critical value, oscillations in shape of the same time scale occur 
before the steady state is achieved. 
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